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pi I Magnetic properties of quark matter are discussed in the hght of the observation of pulsars. 

^^ i Our works about spontaneous spin polarization and spin density wave are reviewed and their 

f^ I implications on compact-star phenomena are discussed. In particular, the former subject may 

^ii: , be directly related to the origin of strong magnetic fields. An inhomogeneous state emerges 

r I I following the chiral transition, where a kind of spin density wave develops. 
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10.1 Introduction 



c^ 



l^ • Pulsars are rotating neutron stars and have a strong magnetic field, by which they emit optical, 

^ • radio , X-ray and 7-ray pulses. 

Although it is not well understood yet about where or how the pulsar emission takes place, 
the magnetic dipole model of pulsars simply states, from the viewpoint of energetics, how the 
rotation energy of neutron stars is converted into pulsar emission [ilj. First of all, let us briefly see 
the physical mechanism. Consider the sphere with radius a composed of uniform magnetization 
M (M//z). Then the solution of the vector potential is 

A^ = —Ma'^—sm9, (10.1) 

3 r> 

where (r>, r<) are the larger and smaller of (r, R) [2\. The magnetic field B is then given by 

1 d 

r sm 6 od 
1 d 

r or 
B^ = 0. (10.2) 
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The magnetic field outside the sphere is 
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Hence it takes the maximum magnitude at 9 = 0, 

51°"*) = 2|m|r-3, 



:io.3) 



(10.4) 



with the total magnetic moment m = 47r/3a'^M. When we consider a star with radius R {a = R), 
|m| can be related to the magnetic field at the magnetic pole of the star, Bp = 2\in\R^^, by 
way of Eq. (jl0.4p . Accordingly the magnetic- field lines can be drawn by solving the differential 
equation, 
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to give 

outside the star (Fig. I1U.1|) . 



r = R{sm.9 / svhOq) , 



(10.6) 





Figure 10.1: Dipole magnetic field- 
lines outside the star. 



Figure 10.2: Angular distribution of the luminosity. 



The magnetic field is uniform inside the star, 

_B(m) ^ 



2m 



:io.7) 



The light cylinder is defined as a cylinder with an axis along the pulsar rotation axis and 
with the radius re = c/uj, where the velocity of the co-rotating frame with the pulsar attains the 
speed of light. Then the emission activity of the pulsar has been considered to originate from 
the magnetosphere with the radius r^., while its mechanism is not clear yet [T]. 

Denoting the inclined angle of the magnetic moment m to the rotation axis es by a, we 
can evaluate the luminosity (energy radiation per unit time) Lm- Then we can immediately 
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see that the component perpendicular to e^, m±, is responsible to the luminosity, which is 
rotating in the plane perpendicular to the rotation axis, m_|_ = |msina|[ei cos{ujt) +e2sin(a;t)], 
with ej • Bj = 6ij. Thus the luminosity becomes equivalent with the one brought about by two 
oscillating magnetic dipoles (see Fig. 110.2]) . A standard formula gives 

for differential luminosity, where O^j is an angle relative to the rotation axis [2]. 
Finally we have 



2w^ , ,2 

= — ^-rnsma 

B^R^uj^ sin^ a 

= ^^^ — ■ ('"■'« 

for total luminosity. It would be interesting here to see that the luminosity (jl0.9p can be recast 
into the form, 

fBr,R^V 9 S(rc)2 , 

Lrn^c[^^] r2 ~ c-^ ■ 47rr2, (10.10) 



. 2r3 y '^ Svr 

which implies that the energy density of the magnetic field (~ B{rcf' /Sir) fiows outward through 
the surface of 47rr^ with the velocity of light c. 

Assuming that this energy release originates from the rotational energy of neutron star, 

^ = -Lm. (10.11) 

where the rotation energy is given as E^ot = \l^'^ with the moment of inertia I, assuming the 
rigid body. From Eqs, (jl0.9p . (jlO.lip . we have 

2TT'^R^sm^a n 
PP= 3j,3 Bl (10.12) 

where P = 27r/a; is the period of the rotation. The characteristic age Tc is defined by 

Tc = P/2P, (10.13) 

which gives a rough estimate of the age of pulsars. 

The spin-down rate co is given as ci; oc — cj'^ from Eqs. (jlO.Op and (jlO.lip . so that the braking 
index n = —ujuj/{ijj)'^ is three in this model. Note that the observed values of the braking index 
are less than three, e.g. n = 2.515 it 0.005 for the Crab, and it is usually considered that 
the significant fraction of the rotation energy is carried away by a pulsar wind [Ij. Thus the 
magnetic dipole model is not sufficient to understand the emission processes in pulsars, while it 
can provide plausible estimates for the magnetic field and the characteristic ages. 
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Nowadays observation oi P — P for pulsars are summarized on the plane 0] , where we can 
see three clusters of pulsars: ordinary radio pulsars, the millisecond pulsars and magnetars. We 
can easily see that most of the radio pulsars are centered around lO^^G. Magnetars are compact 
stars with huge magnetic field {B = O(10^^~^^)G), including soft 7-ray repeaters (SGRs) or 
anomalous X-ray pulsars (AXPs)[3l [5l[6]. It exceeds the critical magnetic field defined as 

Bqed = ^ ^ 4.4 X lO^^G^ (10.14) 

en 

which is obtained by equating the spacing of the Landau levels, 5Ee = fieB/{mec), with the 
mass, rrieCp. So it presents a criterion about how the relativistic effects are important and the 
quantum effects such as e^e~ creation or photon splitting becomes significant. It should be also 
interesting to recall the relation, IMeV^ ~ lO^^G in the natural unit. So, B = lO^^G expected 
in magnetars may correspond to lOMeV, which gets to the strong-interaction energy scale or 
the Fermi energy of nuclear matter or quark matter at nuclear density. 

We can easily asses how the magnetic field affects EOS inside compact stars. As shown in 
Eq. (110. 7p we can consider the uniform magnetic field inside the stars. Then 6Ep is small for 
protons, 6Ep <^ rup, while it is comparable with light quark mass, 6Eq ~ niq, for B = O(10^^)G. 
On the other hand, the magnetic-interaction energy can be estimated by a simple formula, 
-E-mag = fJ-iB, with the Dirac magnetic moment, /ij = ei/{2mi), for i-th particle having mass rrii 
and electric charge ej. For electrons it gives O(keV) for the canonical value of B = O(10^^)G, 
which is comparable with the atomic-energy scale. So we may easily expect that thermodynamic 
quantities or EOS of the neutron-star envelope or the crust region should be much affected by 
the magnetic field [?]■ On the other hand, it becomes very tiny for protons (-Emag ~ keV) even 
for B = O(10^^)G, which implies that the magnetic field little affect the EOS of nuclear matter. 
Since it amounts to -Emag ~ MeV for light quarks, the magnetic field looks to modify the EOS 
of quark matter. However, the Fermi energy should be very large (O(lOO)MeV ^ SEq), and a 
large number of the Landau levels are occupied, which may be well approximated by the usual 
Fermi sphere. Thus we can conclude the effect of the magnetic field little modifies the EOS in 
the core region, even for B ~ lO^^G. 

Other high- field radio emitting pulsars with B > 10^^ G or rotation-powered pulsar/magnetar 
transition objects have been also observed [8]. These observations may give a hint about the 
relationship between magnetic activity and neutron star spin-down. 

The origin of such strong magnetic field has been a basic but a long-standing problem 
since the first discovery of pulsars, while many people believe that it originates from dynamo 
scenario due to the charged current or inheritance of the magnetic field from the progenitor main- 
sequence stars (fossil-field hypothesis) [9]. This problem becomes a current issue, stimulated 
by the discoveries of magnetars. Due to the dynamo scenario, magnetic field is produced or 
maintained by the rotation or convection of charged fiuid, but seems irrelevant for generation of 
high magnetic field in magnetars. It requires that magnetars be born with very short rotation 
periods of the order of 1-2 ms, and may not be supported by observations by two reasons [6]. 
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First, The combination of high magnetic field and very rapid rotation is expected to impart 
a high velocity to the neutron star. However, up to now, the observational evidence for large 
spatial velocities in SGRs and AXPs is poor. Secondly, a large fraction of the rotational energy 
of a newly born magnetar, a few lO^^erg, is lost due to the strong magnetic braking. However, 
an estimate of the explosion energy of the remnants containing magnetars yields values close to 
the canonical supernova explosion energy of lO^^erg, implying initial periods longer than 5 ms. 

Fossil-field hypothesis assumes the conservation of magnetic flux during the evolution of 
main-sequence progenitor star to a compact star: B = {Rm/R)'^Bm with Rm,Bm being the 
radius and magnetic field of the progenitor, respectively. It should be interesting to compare the 
radius R with . Taking the sun as a typical main-sequence star, (i?0)Max — several thousands G 
and i?0 ~ 7 X 10^*^ cm. Squeezing the magnetic flux from Rq to R, we have R = Rq{Bq/B)^''^ ~ 
10km for usual pulsars with B = lO^^G, which should be consistent with standard neutron 
stars. However, R < several x lO^^^cm for B ~ 10^^~^^G, which should be compared with the 
Schwartzschild radius, Rs = 2GM/c^ = 3{M/Mq) ~ 3 x lO^cm for M ~ Mq. Thus we can see 
that R should be comparable with Rg and may fall below it in the extreme case, if the fossil-field 
hypothesis is applied to the generation of the magnetic field in magnetars. 

There is another possibility: a microscopic origin due to ferromagnetism or spin polarization 
of hadron matter |10| . Actually Makishima suggested the hadronic origin of the magnetic field 
in binary X-ray pulsars or radio pulsars [11], since the accumulation of the observational data 
shows a peak with a narrow width. Sometimes high-density nuclear matter is expected to 
show a resemblance with '^He system, by the scaling argument |12|. where quantum effects are 
essential and it shows the ^P-type superfluidity in some physical conditions. Correspondingly 
neutron matter has been shown to exhibit the ^^2 superfluidity around the nuclear density |12j . 
It is also well-known that '^He is very close to a ferromagnetic state under the high pressure 
due to the large magnetic susceptibility. Microscopic calculations of nuclear matter have been 
repeatedly performed to find out a possibility of spontaneous spin polarization inside pulsars, 
but the negative results have been reported so far |13] . The magnetic susceptibility of nuclear 
matter increases monotonously with density, so that the ferromagnetic phase is not be expected 
in nuclear matter at any density. At high-density we can also expect quark matter and we may 
be tempted to ask the possibility in quark matter. Nowadays there have been actively studied 
about the possible evidences of hadron-quark transition in astrophysical phenomena as well as 
relativistic heavy- ion collisions [H ll4l[T5] . If such phase transition occurs inside compact stars, 
it may affect thermal and magnetic properties as well as equation of state (EOS). Here we ask 
a microscopic origin of the magnetic field by considering uniform quark matter inside compact 
stars [16]. 

We consider the possibility of spontaneous spin polarization by the use of QCD. First, the 
evaluation of the total energy for spin-polarized quark matter has been done by using the one- 
gluon-exchange (OGE) interaction at zero temperature [Ml [T71 [18] . Subsequently, the magnetic 
susceptibility has been studied within the Landau Fermi-liquid theory [19j, taking into account 
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the screening effect for gluon propagation |2Ut ETj 122] . Finally we present a magnetic phase 
diagram on the density-temperatm'e plane. 

We also intend to address another interesting magnetic aspect of quark matter, which is 
characterized by the spatial modulation of the magnetic moment [23j . It is an inhomogeneous 
phase accompanying the chiral transition. In the standard scenario, the quark condensate, which 
is the order parameter of spontaneous symmetry breaking (SSB) of chiral symmetry, decreases 
due to the Pauli principle as density is increased. It is eventually vanished at some density- 
temperature point, which means the restoration of chiral symmetry. However, this may not be a 
unique scenario. Recall the FFLO state in the superconductivity; it has a spatially non-uniform 
order parameter and is considered to appear in the vicinity of the critical point, when two 
Fermi spheres with different spins have different sizes [21] . Recent studies have shown that there 
possibly appear various kinds of inhomogeneous phases in the vicinity of the chiral-symmetry 
restoration |25t \26[ [27] . Among them we consider a special one called dual-chiral-density-wave 
(DCDW) state in detail, because it exhibits many interesting theoretical features and is expected 
to bear various implications in the light of compact-star phenomena. 

Nowadays there have been done many theoretical studies about the deconfinement transition 
in high-density nuclear matter by using the MIT bag model [28] or other effective models of 
QCD to find the EOS including quark degrees of freedom. Unfortunately the lattice simulation 
is not possible, because the numerical calculation suffers from the so called sign problem [29] . 
Thus there are ambiguities about the properties and the critical density of the deconfinement 
transition. Here we only assume the presence of quark matter without resource to EOS and the 
details of the deconfinement transition. 

10.2 Spontaneous spin polarization in quark matter 

A simple idea about spontaneous spin polarization owes to Bloch, who first discussed the emer- 
gence of ferromagnetism in electron gas at low density |30[ [31] . Consider the electron gas in the 
positively charged background to compensate the electromagnetic charge of electrons. Then the 
Coulomb interaction is classically vanished, but the Fock exchange interaction between electrons 
with the same spin gives an attractive effect due to the Pauli principle; in the non-relativistic 
approximation it reads 

E.. = -vie' (A) n^/3 [(1 +p)'/' + (l-p)'/'] , (10.15) 

where n is the total number density and p is the polarization parameter defined hjp= (ns=+i — 
ns=-i)/n with the number densities, Ug, of electrons with the spin s = ±1. This spin dependence 
may be easily understood by observing that the repulsive Coulomb interaction is effectively 
avoided for electrons with the same spin, since they cannot approach due to the Pauli principle. 
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On the other hand the kinetic energy is easily evaluated as 



Ekii 



h'6 

V vr 
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'_iy\^/s 
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(l+p)5/3 + (l_p) 



5/3 



(10.16) 



for non-relativistic electrons. Note that the density dependence of each term is peculiar. Search- 
ing the minimum of the total energy etotai(p) = (^kin + Ecx)/V with respect to p by fixing the 
total number density n+i + n_i = n, we can observe that electron gas is completely polarized 
for n^/'^ < He = a/(l + 2^^/'^) with a = 5/(67r^)e^m/?i^(97r/2)-^/^. The phase transition is of 
weakly first order in this case. 
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Figure 10.3: Energy density (arbitrary P 

scale) as the function of the polarization 
parameter p = (n-|_i — n_i)/ng. 

Figure 10.4: Energy per particle for quark matter 
with the OGE interaction as a function of the po- 
larization parameter, p. 

This is the result within the Hartree-Fock approximation, but it has been also shown by the 
quantum Monte Carlo simulation that the electron gas is in the ferromagnetic phase at very 
low density |32) . Recently it has been experimentally observed [33]. A lesson we learned here is 
that we need no spin dependent interaction in the original Lagrangian to see spontaneous spin 
polarization, but a symmetry principle plays an important role in this context. 

When we consider quark matter within QCD, we notice that the situation seems to be similar 
to the electron gas; quark matter is color neutral as a whole, so that the exchange term of the 
one-gluon-exchange (OGE) interaction gives a leading-order contribution to the total energy. In 
ref. we have calculated the interaction energy of relativistic quarks in a perturbative way, and 
demonstrated the spontaneous spin polarization around nuclear density, n^ ~ /Oq — 0.16fm^^. 

Assuming that spin-polarized quark matter with density Uq inside compact stars, we can 
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roughly estimate the magnetic field at the surface (see Eq. I1U.4|) ). 

JD 



max — 3 V^y A'q"'? 

for the quark core of the radius r^, where fiq is the quark magnetic moment and /Utv the nuclear 
magneton □ . 

10.2.1 Magnetic susceptibility within the Landau Fermi-liquid theory 

In the recent papers we have studied the magnetic susceptibility within the Landau Fermi-liquid 
theory |20[ [2T] , to get more insight about the phase transition and derive a phase diagram in 
the density-temperature plane. 

Applying a tiny and uniform magnetic field B (along the z-axis), we can study the linear 
response of quark matter by calculating the change of magnetization, M, by the external field. 
Then the magnetic susceptibility xm is defined as xm = dM/dBl^^^^. We can easily see that 
Xm measures the spin-spin correlation in the normal quark matter, or the curvature of the free 
energy at the origin with respect to the order parameter, the magnetization M. The free energy 
is given as a function of the external magnetic field, F{B), from which M is given by 

M.-^. (10.18) 

Defining the Legendre transform of F such that 

G(M) = F + M-B, (10.19) 

G satisfies the reciprocity relation, 

Therefore xm can be written as 

'd'^G{M) 



which is the inverse of the curvature of G at the origin, assuming M = for B = 0. Adding 
— M ■ B to the free energy G(M) we find the fluctuation amplitude at T, 

{M^)\b=o = Txm, (10.22) 

which means the fluctuation amplitude diverges as xm 



^One may wonder about the back-reaction of the generated magnetic field to quarks. The Fermi energy is 
O(100)MeV and the quark mass is rather heavy in this case, so that the Landau orbiting and the magnetic 
interaction are safely discarded. 
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Figure 10.5: Two Fermi spheres with different spins have different sizes to generate the magne- 
tization after applying an external magnetic field. 



Xm can be written in terms of the Landau- Migdal parameters derived from the quark-quark 
interaction, 

where go = i|ki=fc c?^k/4'7r(7D(k) is the effective gyromagnetic ratio. N{T) is the effective 
density of states around the Fermi surface, which is written as 

iV(0) 






(10.24) 



at r = 0, where the Fermi velocity vp is given by the spin-independent Landau-Migdal parameter 
/f , kp / Ep — {Nckp /Stt"^) fl . /° is the spin-dependent Landau-Migdal parameter. The divergence 
of Xm or 1 -|- N{T)f°' = implies the phase transition to the ferromagnetic phase. 
We present here a heuristic argument by the use of the naive OGE interaction. 



/:qcd = -W'-G^, + ^ gj (i^^D^ - ruf) qj, 



(10.25) 



where 



D^qf 



d^Ay - duA^ - ig[A^, Ay] 
{df, - igAf,)qf 



(10.26) 



with the gluon field, A^ = X]a=i ^u'^'*/2- Since spin is coupled with motion in relativistic 
theories, we must define spin polarization in a proper way. A relevant one is to introduce 
the space-like four vector a'^ with the constraints, a ■ k = and a^ = —1 for a particle with 
momentum k. The standard choice may be 



a = C + 



m{Ek + m) 



m 



(10.27) 



where the three vector ^ specifies the direction of spin in the rest frame of each particle [35]. 
Then we have the polarization density matrix p(k, (^), 

1 



p{kX) 



2m 



Ht + m)P{a), 



(10.28) 
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with the projection operator, P{a) = (1 + 75(i)/2. Given the quasiparticle interaction fj^^a-ac'b 
between quarks with momentum fc, spin polarization C, and color a and q, C,' and h, the color 
symmetric interaction is generally given as 

1 v^ mm 

JkCqC = J^l^ JkCa,qC'6 = "^-^MkCqC', (10.29) 

with the invariant matrix element within the OGE interaction, 

Mkcqc = V^tr (A„/2A„/2) M^'^{k, (; q, C')^m.(^ " q)- (10-30) 

c 

Here D^^, is the gluon propagator and M^^{k,C,]q,C,') the interaction tensor, 

M^^^{k,Q-qX') = tr [l^ p{k,CW p{q,C')\ , (10.31) 

Taking the Feynman gauge for the gluon field, the invariant matrix element can be explicitly 
written as 

^fcCQC = a'jf^ [2"^' -k.q-w?a-b\ j^—^. (10.32) 

In the non-relativistic limit, it is reduced to 

Thus a pair with the parallel spin (^ = (^') gives an attractive interaction in the nonrelativistic 
limit. This makes a base of the Bloch mechanism [31]: the Fock exchange interaction give rise 
to an attractive effect for the pair of the parallel spin due to the Pauli principle. Then the 
quasiparticle interaction on the Fermi surface is given as 

1 



4c, 



qc 



2N?-l 



\k\=\q\=kp ^ SN^EJ 



c^F 



2m — Ep + kp cos 9j-~ — m a -b 



-4(1- cos 0^) 



(10.34) 

We can immediately see that it diverges for the collinear momenta, k = q. Substituting the 
explicit formula of the spin vector in Eq. (|10.27p . we can see that it consists of two parts, the 
spin-independent (/^ ) and spin-independent (/^ ) terms; 

/kC,qC' = /k,q + CC7£,q, (10.35) 

from which we can derive the Landau-Migdal parameters by the angle-integral over the Fermi 
surface. For the spin-independent Landau-Migdal parameter, ff, 

g^(Nl-l)m^ [\ m cos^ 

the spin-dependent one, f^, 

la^ [d^k fdn^ .a I g\Nl-l)m(2Ep+m) ^^ 

^ ~ J Att J Att "^k,ql,k|=|q|=fc^- 8^r2£;^ 3^ +-^i/"^- ^'^■'^'> 
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Putting Eqs. (|1U.36|) and (|lU.37p in Eq. (|10.23p . we have the final expression, 

Cfg' 



(XAf/XPauli) 



-1 



llTT'^EphF 



m{2EF + m) 



(10.38) 



Af2-1 



with the Casimir operator per color, Cf = ^n ' where XPauU is the susceptibility of non- 
interacting Fermi gas, XPauU = gofJ-qNckpEp/'i-K'^. . Thus we can observe that infrared diver- 
gences cancel each other in the expression of xm- It is to be noted that xm is proportional to 
the quark mass. Thus, we may say that heavier quarks favor the spontaneous magnetization. 
Taking the non-relativistic limit, m^ kp, and replacing Cf and g by one and e, respectively, we 
can recover the standard formula for electron gas interacting with the Coulomb potential|31j. 



(x/xo) 
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em 



(10.39) 



In Fig. I10.2.1l we present a result by using the MIT bag model parameters. 
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Figure 10.6: Spin susceptibility xm at T = within the Landau's Fermi liquid theory. XPauii 
is the Pauli paramagnetism. It diverges around kp — 1.4fm~'^, the order of nuclear density, po, 
and quark matter is in the ferromagnetic phase at lower density side. 



10.2.2 Magnetic phase diagram 

We have also studied the screening effect for OGE interaction on xm, since it is well known that 
such many-body effect is important in the electron gas. Actually the screening effect disfavors 
the spin alignment [36]. Moreover, it is necessary to include it when we consider the gauge 
interaction like QED/QCD, since the infrared behaviour is singular and it is improved by the 
screening effect. In QCD it has another interesting effect, depending on the number of flavors: 
it may favor the ferromagnetic transition for a large number of flavors. In the following let us 
briefly see the screening effect. 

The gluon propagator D^y{p) is generally written as 



D,,{p) = P'Dt{p) + P'Diip)-^ 



P/iPu 
P^ 



(10.40) 
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where P^L is the projection operator onto the transverse (longitudinal) mode, 

Pj^u = (1 - 5'mo)(1 - 9,^0) ( -ff/^z. - -f-wj 

P',u = -9,u + ^-Pl.. (10.41) 

Correspondingly, Di(i^{p) are the propagator of transverse (longitudinal) gluons, modified by 
the medium effect, which can be described by the Debye screening [57]. The self-energies for the 
transverse and longitudinal gluons are calculated by the hard dense loop (HDL) resummation 
to give 

nt(po,p) = -* >^ -. |— T' (10.42) 

— 4 p 

f=u,d,s 

in the limit Po/|p| — ^ 0, with the Fermi velocity upj = tpj/Epj and the Debye mass, m^ r = 
g'^Ufkpj/'^TT'^ for each flavor / |37j |f| . The appearance of the imaginary part in lit physically 
means the damping of the transverse gluons due to the interaction with the surrounding quarks 
(Landau damping). Thus the longitudinal gluons are statically screened by the Debye mass, 
while the transverse gluons are dynamically screened due to the Landau damping. Accordingly, 
the screening effect for the transverse gluons is ineffective at T = 0, where soft gluons (po/|p| ^ 
0) contribute. At finite temperature, gluons with po ~ 0{T) can contribute due to the diffuseness 
of the Fermi surface and the transverse gluons are effectively screened, which gives rise to another 
interesting consequence for XM- We have seen a non-Fermi liquid behaviour, which is inherent 
in the gauge theories: there appears T^lnT term in xm, besides the usual T^ term [20|[2T]. 

Taking into account the screening effect, we can evaluate the magnetic susceptibility. To 
summarize, we present in Fig. 110.81 the phase diagram of ferromagnetic quark matter on the 
temperature-density plane. 

10.2.3 Spin wave 

Ferromagnetic quark matter is specified by the non-vanishing magnetization M, so that rotation 
symmetry is violated there, S0{3) — )• 0{2). Such SSB should accompany the Nambu-Goldstone 
mode, spin wave. Different from the usual treatment of the spin wave in the Heisenberg model 
[38], we must carefully formulate it since quarks freely move there. For the ferromagnetism of 
itinerant electrons. Herring took an intuitive but correct approach called the spiral approach 
|3H 139] . He introduced the collective variables as variants from the spin direction. In ref.[40j 
we have proposed a formulation within the path integral of the spin coherent-state. Since the 
spin wave or its quantized magnons plays important roles for the thermal evolution as well as 



^The Debye mass is given as e^ fj^up/i^^ for electron gas in QED. 
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Figure 10.7: Magnetic phase diagram in the density-temperature plane. The open (filled) circle 
indicates the Curie temperature at kp = 1.1(1.6) fm^^ while the squares show those without 
the T^ In T term. 



the magnetic evolution of compact stars, more elaborate studies are needed about its properties 
including the dipersion relation: specific heat and thermal conductivity may be affected by 
the magnons and a new cooling process may open through the absorption and emission of 
magnons in the interaction vertices. Moreover, the magnon-exchange interaction may work in 
the ferromagnetic phase, which gives, e.g., an exotic superconductivity [41j. 



10.3 Chiral transition and Inhomogeneous phases 

Here we consider another type of magnetism in quark matter. Recently there are many works 
about the inhomogeneous phases accompanying the chiral transition \lb\ [26} 127] . According 
to the usual discussion about the chiral transition, the scalar condensate {qq) is decreased as 
density or temperature is increased, and eventually vanished at some density [H]. Since the 
dynamical mass is proportional to the condensate, it also decreases. 

We, hereafter, consider two-flavor (n, d) quark matter for simplicity. Chiral symmetry is 
then represented by SU{2)l x SU{2)ji algebra, under which the quark field is transformed s.t. 






(10.43) 



for small parameters, a",/3" <C 1. The first one is the usual isospin transformation, while the 
second one, called chiral transformation, mixes the different parity states. Many studies have 
shown that Lagrangian should be chirally invariant for massless quarks, but it is spontaneously 
broken in the vacuum, which is specified by the non-vanishing qq condensate, (V'V') 7^ 0. As 
important consequences, the quark mass can be dynamically generated by the spontaneous 
symmetry breaking (SSB) and pions emerge as a collective state of quarks and antiquarks. At 
finite density the Fermi sea prevent the formation of the condensate due to the Pauli principle, 
and chiral symmetry should be eventually restored at some density. 
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Recent studies have suggested that the homogeneous ground state becomes unstable for 
producing the inhomogeneous phases in prior to the chiral transition. Here we consider one 
of such instabihties in detail by using the effective model of QCD, Nambu-Jona Lasinio (NJL) 
model |231 I42j . In the paper [23j we have demonstrated formation of a density wave in quark 
matter in the chiral limit: where not only scalar density but also pseudoscalar density takes a 
non- vanishing value, 



{■^p^p) = A cos{q ■ r) 
ii'ilbTsi') = Asin(q-r). 



(10.44) 



Note that we can utilize the pseudoscalar degree of freedom in quark matter, different from the 
vacuum, which must be an eigenstate of parity. Consequently the above configuration means 
a violation of parity. It should be interesting here to see a similarity with the FFLO state in 
superconductivity [2jj: a combination of both densities can be written as {tpil^) + i{'4>i'y5T3ip) = 
Aexp(ig-r), which corresponds to the complex energy-gap parameter due to the condensation of 
the Cooper pairs [43j. In the FFLO state, the gap function, A(r) = {ip^{r)tp^{r)), modulates in 
space with the wave number characterized by the center-of-mass momentum of the Cooper pairs. 
The similarity indicates a general mapping between superconductivity and magnetism [44j. It 
might be also interesting to see the phase transition in terms of the susceptibility. Similarly to the 
magnetic susceptibility in the previous section, we can introduce the susceptibility or the density 
correlation function in scalar or pseudoscalar channel, Xs(ps)(<^i q) [23]. In the ferromagnetic case, 
Xm(w = 0,q = 0) diverges at the critical point, but Xs{ps)('^ = 0, g 7^ 0) in this case. The finite 
value of q is characteristic in the inhomogeneous phase. 



<vN^ 




< viy^T-'v) 




Figure 10.8: Dual chiral density wave 
in the chiral space. 



Figure 10.9: Possible chiral restoration paths on the 
A — q plane. Dotted lines show the first order phase 
transitions. 



Consider the 2-flavor Nambu-Jona Lasinio model as an effective model of QCD at low density, 

>Cnjl = Hi^ - m,)i; + G [{^i^f + (V^iTsrV')'] , (10.45) 

where G is the coupling constant and nic the current mass. □ Putting Eq. (|10.44p in the 



*We, hereafter, consider the chiral limit where vac = in this subsection. 
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Lagrangian under the mean-field approximation, we have 



Cmf = tplic^ - M* (cos(q • r) + ^75X3 sin(qr • r))] 



M 



*2 



AG ' 



(10.46) 



where the effective mass M* is defined by way of M* = —2GA. Introducing a new spinor by 
using a local chiral transformation, 



Cmf is recast into 



V;(r) = e-*^"^9-^/V«,W, 



^MF = ipw [i^ - M* - r3757 • q/2] ip^ 



M 



*2 



AG 



(10.47) 



(10.48) 



This is a kind of the Weinberg transformation ^45] , which makes the Lagrangian in the presence 
of inhomogeneous DCDW to the one with the homogeneous axial-vector field. 

Table 10.1: Separation of amplitude and phase degrees of freedom by way of the Weinberg 
transformation. 



non-uniform 



uniform 



{ipip) = A cos q ■ r 
{ipi'y^T^iJj) = A sing • r 



{'4'w'4'w) = A 

(ipwilbTsipw) = 

75T37Q/2 



Then the quark wave function is still given by the plane- wave solution V'fc(^) with the energy 
eigenvalue, 

E^{p) = yjp"^ + M*2 + q2/4 ± y^{p ■ qf + M*^. (10.49) 

Accordingly the Fermi seas are deformed in oblate and prolate shapes, depending on the spin 
degree of freedom. In Fig. we show energy spectra for example. 

Thermodynamic potential is given as a function of M* and q at finite temperature T and 
chemical potential /j, : 



n{q,M*) 



-TlogZ/V 



-NfN, 



(2^ 



s=± 



-p{E^-^l) 



+ 1 



-/3(£;-+m) 



+ 1 



+ E'\ + 



*2 



^.50) 



The last term in the curly bracket denotes the vacuum contribution, which has an essential 
role in the chiral transition. Since single-particle spectra are deformed in momentum space, the 
simple momentum cut-off is irrelevant and sometimes leads to unphysical results. We employ 
the proper-time regularization with cutoff A to evaluate the vacuum contribution. 
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Figure 10.10: Energy spectra for p± = under the condition, q/2 > M*. E is the one in the 
absence of DCDW. 
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Figure 10.11: Phase diagram of DCDW on the density-temperature plane. DCDW emerges in 
the vicinity of the chiral transition. 



10.3.1 "Nesting" mechanism 

First we consider the one dimensional case. Assume the presence of the spin density wave (SDW) 
exp(ibigz) coupled with the non-relativistic electrons like 



Then eigenvalues are 



Hu 



^k = -^(e/c + f-k+q) ± 



y(e-^''^uV+u + e^''^nV- 



(efc - eu+qf + W 



1/2 



(10.51) 



(10.52) 



with the free particle energy, e^ = k"^ /2m (Fig. ll0.12|) . In the case of g = 2kp, there is produced 
a gap at the Fermi surface of the free electrons. Thereby the total energy is always decreased 
by the interaction, independent of the strength of the coupling. This is called a nesting effect 
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of the Fermi surface [461 W7\ . We can also see this by calculating the correlation function; it is 
given by the Lindhard function L{uj = 0, q), which exhibits a logarithmic divergence at g = 2kF- 
In the higher dimensional cases, the nesting is incomplete but we can see its reminiscence. 




kF q=2kF 

Figure 10.12: Energy spectra of one dimensional electron gas in the presence of SDW. 



We can see the similar mechanism should be also responsible to the emergence of DCDW, 
but in somewhat different manner. By putting p_|_ = 0, we can consider quasi-one dimensional 
case. The energy spectra exhibit different features, depending on the conditions, q/2 > M* and 
q/2 < M* . It is not obvious which condition holds in the DCDW phase, but the numerical 
results shows that the former condition always holds. Then we can understand the mechanism 
of the DCDW formation by two steps. First looking at Fig. 110. 12|, we can see that the levels of 
massless quarks cross each other at pz = 0. After switching on the mass term the level crossing is 
finally disentangled to give E^. So, if we set q = 2fj,, there is always the energy gain, irrespective 
of the dynamics. 

10.3.2 Some implications 

The presence of DCDW has some implications on compact star phenomena. Let's consider, for 
example, the /3-decay of quarks in the DCDW state, 



d{pi) -^ u{p2) + e (ps) + Ve{pi). 
From Eq. (jl0.47p the DCDW state can be represented as a chirally rotated state, 

IV'^) = ^(q)|V), 

with 

U{q) = exp(i / A[]q • rcf'x), 

the quark current /ij'+j2 = '07'^(1 — 75)t_|_'0 is transformed as 



h'^ 



l+«2 



^i^l) K+i2^(.l) = exp(iq • r)/i^_ 



1+J2' 



(10.53) 

(10.54) 
(10.55) 

(10.56) 
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Thus we can see that DCDW supply an extra momentum to modify the momentum conservation 
in the /3-decay process. It is easily shown that the neutrino emission due to the quark /3-decay is 
suppressed by the energy-momentum conservation: the momenta of all the particles should be 
collinear for reaction (J10.53P for free massless quarks at low temperature. So DCDW catalyzes 
the reaction (jl0.53p . which may give rise to a fast cooling of compact stars [48J. 

The symmetry breaking pattern in the DCDW phase is as follows: the original translational 
symmetry along z-axis (the generator is pz) and U{1) subgroup of chiral symmetry (the generator 



is Q^) is broken at the same time, but it is still invariant under a combination of these operations, 



^P. Xt^Q3(l)^^p.+Q3(l). (10.57) 

Actually we can see that the single-particle wave function is the eigenfunction of the operator, 

exp{ipza) exp(-i75r3/26')'0fc(r') = exp{ikza)tpk{r) (10.58) 

with 6 = qz. The collective modes ("phasons") have then hybrid properties of phonons and 
" pions" . We can derive the dispersion relation of the collective modes by way of the Ginzburg- 
Landau theory. If we consider the phase fluctuations, u(r, i), the free energy should be written 
as 

F = Fo{A,q)+6F{u{r,t)), (10.59) 



with 



,..1 



A{Vzuf + B{Vluf + c(^^' 



d^r, (10.60) 



like liquid crystal (smectics) [12], where the coefficients A,B,C are the function of A, q^\- The 
Euler-Poisson equation for u then reads, 

C^-A^ + B{Vlfu = 0. (10.61) 

The plane-wave solution, u = uqc*''^*^ \ is obtained with the anisotropic dispersion relation, 

w^ = {Ak'i + Bki)/C. (10.62) 

Note that the dispersion relation exhibits a hybrid nature of type I and II Nambu-Goldstone 
bosons [50j: oj (x kz for the longitudinal excitation with k± = 0, while w oc /c^ for the transverse 
excitation with kz = 0. 

10.3.3 Magnetic aspect of DCDW 

Here we reveal the magnetic aspect of DCDW [23]. Since explicit form of the spinor can be 
analytically obtained within the mean-field approximation, one can evaluate various expectation 
values in the DCDW phase. 
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By the use of the spinor ipy^ one can evaluate the expectation values of bilinear form of 
gamma matrices O in a simpler way, 

Spin expectation value, for instance, O = 707573/2 = Sz, vanishes because it is proportional to 
the stationary condition with respect to wave number q, 



■{i'l{p)^z1pw{p)) OC 



(27r)3^^"'^^^^^"'^^^' dq 



= 0. (10.64) 

DCDW 



On the other hand, there is formed a spatial modulation of magnetic moment along with 
DCDW. The Gordon decomposition of the gauge coupling term provides the magnetic interaction 
in presence of external gauge field: ^^pr'tp(7fj_uipF'^'^. Only z component of the anomalous magnetic 
moment remains finite, 

mr)a23Hr)) = {i^{r)a3Mr)) = 0, (10.65) 



{^p{r)al2'lp{r)) = / 77^:^3 (V'l(p)7oo-i2^«,(p)) cos(g • r) 



where 

2M* 
{'^i{p)lQ'yi2^w{p)) = ,^^„ „ [n+(p) - n_(p)] (10.66) 

with n± being the Fermi-Dirac distribution function with energy spectra e±. Therefore, the 
magnetic expectation value reflects the anisotropy in momentum distribution of the DCDW 
phase, and depends on finiteness of the effective mass M* . 

To see a possible consequence of the spatially modulated magnetic order, we evaluate its 
local magnetic-flux density with effective charge Q = (3 — i) e, 

Q (Vi(r)ai2^(r)) ^ 
''=2AF (^t(,)^(,)) 3^^' (10-67) 

where we have evaluated magnetic expectation value par quark, and multiplied by quark density. 
For densities normalized by the normal nuclear density Pb/ P ~ 3—4, the local flux $ is estimated 
to be 0(10^^) G, which well exceeds the critical magnetic field -Bqed- 

It might be interesting to recall that the stable Hartree-Fock solution for the case of the short- 
range potential is always either the normal paramagnetic state or the uniform ferromagnetic state 
|31j . We have considered the zero-range interaction between quarks within the NJL model, but 
spin density wave emerges in our case. 

10.3.4 Deformed DCDW 

We must take into account the symmetry breaking effect or the finite current mass for the 
extension of the above framework [51]. It is important for a realistic discussion of the phase 
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transition at moderate densities. Moreover, it may become important to elucidate the appear- 
ance of various inhomogeneous phases. The real kink crystal (RKC) can be easily generalized 
to include massive quarks, and has been studied on the temperature-density plane [26]. Their 
results show an interesting change of the domain of the RKC as the current mass is changed; it 
looks to shrink for larger mass. 

Recently we have shown that the symmetry breaking effect can be taken into account by 
a variational approach, without spoiling the original features [T0|. Using the ansatz for the 
condensates, 

{'ipil}) = Acos0(r) 
{i'ilbni') = A sin 6'(r), (10.68) 

represented by the spatially dependent chiral angle 0{r), one must solve the Hartree equation 
self-consistently in the presence of current-mass term, but it should be a hard task. Instead, we 
apply approximate methods, like a perturbative method [52]. Here we take a variational method, 
which may give better results than the perturbative one. Anyway the symmetry breaking effects 
should be small in the light of the success of the chiral symmetry approach to various phenomena. 
So we take into account only the deformation of the chiral angle, leaving the quark wave function 
unchanged. Then we can see that must satisfy the sine-Gordon equation, 

(fOiz) 
Ca ,\ -FmcAsin6'(z) = (10.69) 

at the leading order of rric. It should be interesting to observe the appearance of the sine-Gordon 
equation. The relevant solution is then written as 

e = T: + 2axn{mlz/k,k), (10.70) 

in terms of the Jacobian elliptic function, where we introduced the effective pion mass by the 
relation, Catti^ = — w-c^; and am(^, k) is the amplitude with modulus k. In the limit, A; ^ 
with keeping m*^/k =const.(= g/2), — > vr -|- gz so that we have the DCDW solution again in 
the chiral limit. On the other hand, 9 — )• 4tan~^ \e'^^^^ , in the opposite limit. A; — >■ 1, which is 
the well-known kink solution. Thus our variational ansatz may be regarded as an embedding 
the sine-Gordon kink crystal in 1 + 1 dimension into 1-1-3 dimensional quark matter; actually 
we can easily see that Ca — > /^ {f-n'- pion decay constant) and Eq. p0.69p is reduced to the 
standard sine-Gordon equation in 1 -|- 1 dimension, 

^^-m2sin0(z) = O, (10.71) 

by way of the Gell-Mann-Oakes-Renner relation, /^m^ = —rric (vac|gg|vac) [53] . 
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10.4 Summary and Concluding remarks 

It should be important to study EOS at high-density region from various viewpoints. Here we 
tried to extract information about the properties of hadron or quark matter and EOS, con- 
fronting the magnetic phenomena of pulsars. We have seen two kinds of the magnetic properties 
in quark matter: one is spontaneous spin polarization and the other is the spin density wave. 
The former is similar to the itinerant electrons in QED, while the latter is related to the chiral 
symmetry of QCD. These discussions may be in a rather primitive level and must be verified by 
observations or more elaborate theoretical studies toward more realistic description. Compre- 
hensive study of these magnetic properties is also needed to study their competition. 

Energetics or the mechanism of the giant flares observed in some magnetars may give a 
direct hint about the origin of the magnetic field. Thermal evolution of compact stars should be 
important for us to get the information about the EOS and properties of high-density matter. 
Actually the magnetic properties is not directly related to EOS but may manifest through the 
thermal activities of compact stars. 

There have been done many works about color superconductivity (CSC) in high-density 
matter [54J. It also has implications on compact star phenomena. So it should be interesting 
to study the interplay of CSC and magnetic properties in quark matter. In ref.[55j we have 
discussed a coexistence of CSC and ferromagnetic order. 

More elaborate studies are necessary for inhomogeneous phases: relations among the var- 
ious phase should be figured out as well as their properties. For example, there seem many 
resemblances between pion condensation in hadron matter [56] and DCDW state. One may find 
a kind of hadron-quark continuity across the deconfinement transition. They may also appear 
during relativistic heavy-ion collisions. So it is interesting to consider how we can observe them 
in this context. 
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